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Abstract 

We study the zero temperature static properties of dissipative 
ensembles of quantum Ising spins arranged on periodic one dimen- 
sional finite clusters and on an infinite chain. The spins interact 
ferro-magnetically with nearest-neighbour pure and random couplings. 
They are subject to a transverse field and coupled to an Ohmic bath of 
quantum harmonic oscillators. We analyze the coupled system using 
Monte Carlo simulations of the classical two-dimensional counterpart 
model. The coupling to the bath enhances the extent of the ordered 
phase, as found in mean-field spin-glasses. In the case of finite clus- 
ters we show that a generalization of the Caldeira-Leggett localization 
transition exists. In the case of the infinite random chain we study 
the effect of dissipation on the transition and the Griffiths phase. 
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1 Introduction 



The influence of quencfied disorder ^ |2] on tlie one liand and tlie effect of 
tfie coupling to an environment jHl E], on the other, on finite dimensional 
quantum spin models have been separately analyzed recently. The former 
introduce Griffiths-Mc Coy singularities that are particularly important for 
quantum systems PflEl lEj • Dissipation implies decoherence and the coupling 
to a quantum bath generates highly non-trivial localization phenomena at 
least for a single two- level system |7j. The combined effect of quenched 
disorder and the coupling to a quantum reservoir has only been analyzed in 
mean-field models so far [Hlin]- 

A ring of quantum Ising spins is probably the simplest finite dimensional 
model exhibiting a phase transition from a paramagnetic (pm) to a ferro- 
magnetic (fm) phase where one can analyze the mixed effect of disorder 
and dissipation. Indeed, the bath introduces long-range ferromagnetic in- 
teractions in the imaginary-time direction that add to the nearest-neighbor 
ferromagnetic couplings due the quantum nature of the spins and, at zero 
temperature, one finds a phase transition even for finite size clusters. The 
d + 1 equivalent classical problem is characterized by perfect correlation in 
the imaginary-time direction and the effect of disorder is consequently very 
strong. For sufficiently large system sizes and, in particular, in the thermo- 
dynamic limit rare regions in the sample may have a very strong effects on 
the properties close to the phase transition, giving rise to very interesting 
Griffiths phenomena. 

The zero-temperature critical behavior of the isolated finite dimensional 
infinite system in one, two and three dimensions has been studied in detail 
in a series of analytic [H [101 E] and numerical papers [121 113 HH UHl HEl HZ] • 
An 'infinite randomness fixed point scenario', in which the system appears 
more and more random at larger and larger scales, emerged pQ. The critical 
behavior is quite uncommon, with activated scaling on the transition and 
very large distributions of local properties in a finite interval close to the 
quantum critical point that lead to different typical and average values and 
the divergence of the global susceptibility well within the paramagnetic phase. 

In this paper we investigate the effect of Ohmic dissipation on the equi- 
librium properties of small clusters and an infinite chain of quantum Ising 
spins with ferromagnetic exchange interactions of pure and random type. 
We focus on the following questions: (i) How is the Caldeira-Leggett local- 
ization transition modified when several ferro-magnetically interacting spins 
are coupled to an environment, (ii) Which are the characteristics of the quan- 
tum phase transition, as a function of the transverse field and the coupling 
strength to the bath, for finite clusters and in the thermodynamic limit, (iii) 
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Whether the pecuhar phenomenology of the isolated random Ising chain 
especially in its paramagnetic phase, also exists when it is coupled to an en- 
vironment. Conflicting predictions, based on phenomenological droplet-like 
arguments [IHj and the analysis of some simplified models [TH] , appeared in 
the literature recently. We address these questions here using Monte Carlo 
simulations of the equivalent classical model in d = 2 with short-range inter- 
action in the spatial direction and long-range interactions in the imaginary- 
time direction as introduced by the coupling to the bath. 

The paper is organized as follows. In Sect. El we define the model. Sec- 
tion 01 is devoted to a brief description of the Monte Carlo method used to 
generate the numerical data. We also summarize in this Section the parame- 
ters studied. In Sect. HI we explain our results. Finally, in Sect. El we present 
our conclusions and discuss some areas of applicability of our results. 

2 The model 

The Hamiltonian of the random quantum Ising spin chain is 

TV N N 

Hj = -Y: J^m+l - E r.^f - E . (1) 

i=l i=l i=l 

The spins lie on the vertices of a periodic one dimensional lattice with N 
sites. They are represented by Pauli matrices satisfying the SU(2) algebra: 
[o"j'^,(3"^] = SijeabciO'i with a,b,c = 1,2,3 associated to x,y,z, respectively 
and i = 1, ■ ■ ■ , N. We consider Jj, the strength of the exchange interactions 
between nearest-neighbors, as independent quenched random variables in 
the interval < Jj < 1 with uniform probability density. Without loss of 
generality, we choose Jj positive since for one dimensional lattices the sign of 
the interactions can be gauged away (even when the system is linearly coupled 
to the coordinates of the harmonic oscillators representing the bath). The 
next-to-last term is a coupling to a local transverse field that, for simplicity, 
we choose to be the same on all sites, Fj = F. The last term is a coupling to 
a longitudinal field that one may include to compute local susceptibilities. 
The chain plus environment is modeled by 

H = Hj + Hb + Hi + Hct (2) 

where Hb is the Hamiltonian of the bath, Hj represents the interaction 
between the chain and the bath and Hct is a counter-term that serves to 
eliminate an undesired mass renormalization induced by the coupling pUj . 
We assume that each spin in the system is coupled to its own set of N /N 
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independent harmonic oscillators with the total number of them. The 
bosonic Hamiltonian for the isolated reservoir is 

(3) 

The coordinates xi and the momenta pi satisfy canonical commutation rela- 
tions. For simplicity we consider a bilinear coupling, 

N N 

Hi^-Y,^tY.(^ii^i ^ (4) 

1=1 1=1 

that involves only the oscillator coordinates. The counter term reads 

N -1 / N \ 2 

The partition function of the whole system for a particular realization 
of the random exchange interactions is Zj — Tre~^^, involving a sum over 
all states of the chain and the bath. This problem can be mapped onto 
a classical model by using the Trotter-Suzuki formalism that amounts to 
writing the path-integral for the partition function as a sum over spin and 
oscillator variables evaluated on a discrete imaginary-time grid, Tt — tAr 
with At = Ph/Nr, labeled by the index t — 0, • • • , A^^ — 1- We assume 
periodic boundary conditions in all directions with f3h the length in the r 
direction. This mapping is exact in the limit Nt- ^ oo. In what follows, we 
use units such that ks = h = 1. 

The integration over the bath variables can be performed explicitly; for 
an Ohmic bath the resulting classical action reads 

Af^-l N Nr-1 N 

A[K„B,a;s^ = ^ E E - ^ E E ^^^^^ 

t=0 i=l t=0 i=l 



2 ^'h^Nj sin' {n\t-t'\/N^) ' 

where we set hi — and we dropped an irrelevant constant factor, sj are 
N X Nj. classical Ising variables, s* = ±1, representing the z component 
of the quantum spin at each instant at site i. Notice that the bath is 
responsible for the appearance of long range ferromagnetic interactions in 
the T direction controlled by the chain-bath coupling parameter a. From the 
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quantum nature of the spins, a nearest-neighbour ferromagnetic couphng also 
arises in the r direction. Its strength B is defined by 



exp -2B = tanh(Arr) , 



(7) 



and describes the intensity of the quantum fluctuations. Further, since Jj are 
quenched random variables, the random values of the spatial interactions 



are the same along each imaginary-time direction, i.e. depend only on the 
spatial label i. 

3 Monte Carlo simulations 

We analyze the static properties of the system by means of Monte Carlo 
simulations applied to the effective (1 + 1) dimensional partition function. 



The sum over all spin configurations accounts for the statistical average that 
we denote, henceforth, with angular brackets. In the disordered cases, the 
thermodynamic properties and, in particular, the phase transition follow from 
the analysis of the free-energy density averaged over the probability distribu- 
tion of the random exchanges. We indicate this average with square brackets. 

Taking advantage of the fact that the ferromagnetic models are not frus- 
trated (all Ji > 0), we use an efficient Monte Carlo algorithm that adapts 
the cluster updates [21] to the case with long range ferromagnetic interac- 
tions |22]. In addition, in order to avoid running several simulations with 
similar parameters, we use histogram methods to scan the space of parame- 
ters near a phase transition 

The classical counterpart model is defined on a rectangular lattice with 
size N X Nt-. The zero temperature limit, /5 — oo, is achieved by taking the 
thermodynamic limit in the imaginary time direction, Nt- oo. We used 
various system sizes varying the aspect ratio of the 2d sample. We took 
form 1 to 32 and form 4 to 1024. We found that up to 400 Monte Carlo 
sweeps were needed to equilibrate the largest sample sizes (note the efficiency 
of the cluster algorithm!). We simulated up to 2048 different realizations of 
the exchange interactions to account for the effect of disorder. 

We concentrate on the critical properties of finite and infinite periodic 
chains at zero temperature. As we have already set the disorder scale < 



Ki = ArJi 





(9) 



s*=±l 
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Jj < 1, a complete description of the state of the system is given by the 
dimensionless parameters a and ArF. For comparison, we also simulate the 
non-random chain coupled to an Ohmic bath replacing the random variables 
Jj with uniform distribution between and 1 by a constant value equal to 
the maximum, Jj = 1. 



4 Results 

In this Section we present the results of our numerical simulations. 

As mentioned before, we used system sizes ranging form = 1 to = 
32. We separate our results in two groups, small A^ (A^ < 4), and large A^ 
(A^ > 4). In the first group we intend to study whether the random exchange 
interactions modify the physics of the well known Caldeira-Leggett model 
while for the second group our results generalize those reported in ^ (where 
no random interactions are present) and those reported in El ^1 12^ 
(where no bath is present). 



4.1 Phase diagram 

4.1.1 Analysis of the Binder parameter 

The phase transition in the (ArF, a)-plane can be found by analyzing the 
Binder ratio 



9av — 2 



3-7^4 ' (10) 

where the angular brackets indicate the statistical average and the square 
brackets denote the average over disorder (if there is no disorder, this step is 
not necessary). The magnetization density is defined as 

(11) 



i=l t=0 



The analysis of the Binder ratio presents different aspects whether we con- 
sider a finite cluster or its thermodynamic limit. 

(z) Small N 

When A^ is finite the zero-temperature limit of the quantum problem 
corresponds to a classical model defined on a stripe of finite width (A^ < oo) 
and infinite height (A^^ oo). In this case, the usual finite size scaling 
analysis of the adimentional Binder parameter can be applied. The critical 
values are easily found by plotting Qav as a function of a for several values 
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Figure 1: (Color on-line.) The Binder cumulant as a function of a for = 1 
with AtcFc = 0.48. The localization transition occurs at etc = 1.10 ± 0.02. 

of A^^. The dimensionless combination of momenta in Eq. pO|) yields Qav 
independent of size at criticality. Then, ac can be located by finding the 
intersection of the Qav curves against a for several A^,- |25i| for every choice of 
AreFe. 

The case iV = 1 corresponds to the well known problem of a two level sys- 
tem coupled to a bath. In this case we expect a transition from an incoherent 
state to a locahzed one at = 1.0 as analytically demonstrated by Leggett 
et al. 12 . As a test for our numerical algorithm, we re-derive this result 
numerically. In Fig. ^ we show an example with A^ = 1 and AtcFc = 0.48, 
where we see that the curves cross at ac = 1.10±0.02 (the error measures the 
precision on the crossing point). We build the phase diagram by repeating 
these steps for several values of AtcFc- 

We use the same method to obtain the critical line for N = 2 and 4 both 
for random and non-random exchange interactions. In Fig. |2l we show the 
phase diagram for A^ = 1,2,4. In general, it is hard to determine ac for 
ArF < 0.2 where Qav varies less than 0.1% and ac is masked by the noise. 

For A^ = 1 the critical values (AtcFc, ac) are well fitted by a linear func- 
tion that for AtcTc approaches the analytical result Oc = 1 (ZI- The 
localization transition moves to lower values of ac when A^ increases. In ad- 
dition, for finite A^, the disordered chain has a lower AtcFc for a given ac- 
However, this difference tends to reduce when AtcFc 0. 
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N = 4 I — ■ — I 



0.2 0.4 0.6 0.8 1 

AtT 

Figure 2: (Color on-line.) Critical boundary (AtcFc, Oc) for small clusters, 
< 4. For = 1 this corresponds to the localization transition in the 
Caldeira-Leggett model. The data are well fitted by a linear function reaching 
Q^c = 1 at AtcFc = 0. For > 1, we show the transition curve for systems 
with random exchange interactions (thick lines) and non-random exchange 
interactions (thin lines). 



(ii) Thermodynamic limit 

In the thermodynamic limit the finite size scaling analysis must be done 
using the parameters N and A^^ and, naively, one is forced to assume a scaling 
relation between them. The Binder cumulant is an adimentional quantity 
that on the phase transition must scale as 

Qa, ~ giN/^,y) (12) 

where ^ is the spatial correlation length, ^ is a scaling function and y is a 
ratio between the imaginary-time size A'^^ and some adequate function of the 
spatial size A^. For instance, one has 

{Nt-/N^ conventional scaling , 
N^-fe^'' activated scaling . 

A simple argument shows that for fixed A^ and generic values of the other 
parameters the Binder cumulant attains a maximum as a function of A^,- jT2l 
IE] • This argument still holds when long-range interaction are introduced by 
the bath. Indeed, when A^,- is very small with respect to A^, one effectively 
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has a very long one-dimensional stripe (along the spatial dimension) and the 
long-range interactions induced by the bath on the imaginary-time direction 
are irrelevant. Moreover, for the values of the parameters that are close to 
the phase transition of the infinite 2c? system, one is well above the transition 
of the one dimensional system and Qav — 0. In the opposite limit of Nr ^ N 
one is back in the kind of finite cluster problem discussed in the previous 
paragraphs. Now, for the parameters chosen, the new one- dimensional model 
should have a finite correlation length away from its critical line and Qav —>■ 
as well. For and all other parameters fixed Qav should then reach a 
maximum as a function of Nt-. The maximum value g^""^ is independent 
of at criticality, see Fig. El We extrapolate the thermodynamic critical 
values by considering a collection of chains with sizes ranging from = 8 
to = 32. For each At^Tc we looked for the value of a such that the A^- 
independent maximum is reached. The phase diagram in Fig.|3]has been built 
repeating these steps. The error-bars estimate the dispersion on the values 
of Qav arising from the Monte Carlo method and the disordered interactions 
after performing error propagation on Eq. ()10|). 

In Fig m we see that the infinite chain has a phase transition even at 
a = 0, in contrast to the small A^ case. We find that the phase transition 
curve may be accurately described by a power law, and the size of the ordered 
phase increases with increasing coupling to the environment. For comparison 
we also show the transition line for a non-random Ising chain with Ohmic 
dissipation [3] as determined by using the same procedure and similar sizes 
as the ones we use in the disordered problem. 

The critical values for AtcTc with ac = 0, which corresponds to switching 
off the bath, agree with well established analytical results for both the random 
and non-random cases 1211 • Indeed, the critical values of the random 
and non-random models in terms of our variables are AtcFc = 0.25 and 
AtcTc = 0.44 respectively, which agree with our numerical calculations |27j . 



4.1.2 The global linear susceptibility 



In a usual paramagnetic to ferromagnetic phase transition one can expect to 
identify the critical line by looking for the location of the divergence of the 
global linear susceptibility defined as 



N N 
1=1 1=1 



d{mi) 




dhi 


hi=0_ 



1 ^ 



T 1 = 1 



and 



rrii 



E' 

t=i 



(14) 



(15) 
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Figure 3: (Color on-line.) Typical behavior of the Binder cumulant at fixed 
AtcFc and in a neighborhood of its corresponding ac- The maximum of Qav 
is independent of only at a = ac, see the central panel. 
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Arr 

Figure 4: (Color on-line.) Phase diagram for N = oo. We show this line for 
the disordered chain with a uniform distribution of exchanges on the interval 
[0, 1] (thick lines) and the non-random chain with J = 1 (thin lines). 



With the aim of later studying the effect of dissipation on the distribution 
of local susceptibility of the random chain, we attempted to analyze the 
divergence of the global susceptibility in finite clusters |28j. Indeed, when 

< cxD we expect to find a divergence of x only on the phase transition, 
and not within the paramagnetic phase as found for the isolated random 
infinite chain P]. For finite values of the number of imaginary-time slices 
Nr we found that the global linear susceptibility at a fixed has a peak at 
value of ArF that is quite higher than the one obtained from the analysis of 
the Binder parameter in the limit N^- ^ oo. As expected the height of the 
peak increases while its position moves towards lower values of ArF when 
A^^ increases. For A^ = 4, A^,- = 128 and a = 0.45 we find a deviation of the 
order of 10%. This indicates that the finite size effects in A^,- are still very 
important. Note that this holds for the ordered and disordered system as 
well. These observations are illustrated in Figs. El 

In Fig.iniwe attempted to compare the form of the peak for similar ordered 
and disordered systems, i.e. with the same size and under the effect of the 
same external bath. We note that the peak is more pronounced in the ordered 
case and the deviation between its location at finite A^^- and the asymptotic 
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Figure 5: The average global linear susceptibility x for etc = 0.45 as a function 
of AtT. N = 4 and Nr = 128, 64, 32 from top to bottom. Top panel: 
Random interactions. Bottom panel: pure ferromagnet. 

Nj. —>■ oo one is less pronounced. More precisely, we found 

Binder Peak in x 

K-^oo AT, = 128 

0.77 0.84 ordered ^ ' 

0.55 0.66 disordered 

In conclusion, we found that the extrapolation of the finite N^- suscepti- 
bility data to the thermodynamic limit is indeed very hard jSHI- A detailed 
analysis of the scaling behavior needs to be perform but it currently goes 
beyond out computational power. This failure has to be kept in mind when 
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Figure 6: Comparison between the global linear susceptibility for the model 
with random (left) and pure (right) systems with N = 4, N^- = 128 and 
ac = 0.45. 

trying to analyze the dynamic properties of the dissipative disordered Ising 
chain in the paramagnetic phase, see Sect. 14.21 

4.2 Critical scaling: conventional or activated? 

One of peculiarities of the critical behavior of the isolated quantum Ising 
chain with random interactions is that the critical scaling is activated instead 
of of conventional power law type. We wish to investigate whether a similar 
behavior persists when Ohmic dissipation is included. 

4.2.1 The Binder ratio 

The Binder ratio not only gives a criterion to find the phase transition curve 
but also helps to derive scaling laws. The study of the scaling laws for the 
function Qav can in principle provide an answer to the question of whether 
the critical behavior is of conventional or activated type. Finite size scal- 
ing implies that dimensionless quantities should scale as functions of N/C,, 
and as either Nt/N^ in the case of conventional scaling or Nr/e^' in the 
case of activated scaling. At the critical point, the correlation length in the 
spatial direction, ^, diverges. This suggests that the scaling variable should 
be Nt-/N^"-^ or In N^-f In N^"'^ in each case respectively, where iV^"^ is the 
value of Nt that maximizes Qav IIB] ■ This way of analyzing the data has the 
advantage of not needing to determine the critical exponents z or z. 
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7: (Color on-line.) Test of conventional scalings at {AtcTc 
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Figure 8: (Color on-line) Test of activated scalings at {AtcTc = 0.48, etc 
0.25). 



However the interpretation of the numeric results is delicate. Even for 
a = 0, where it is possible to show analytically the activated scaling pp, 
it has been hard to distinguish between both types of scaling by means of 
Monte Carlo simulations IT!^ 1^ . In particular, in ci = 1 the numerical 
data confirmed the analytic prediction only when a mapping to fermions was 
used |T3j. 

In Fig 13 we show the scaling of the Binder parameter, Qav/g^"'^, with 
^l^jnax (conventional) and in Fig El with InA^/lnA^™"^ (activated). Both 
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scalings give almost the same critical values (AtcFc, Oc) in agreement with 
our previous simpler analysis using the independence of the height of the 
maximum on at criticality. Unfortunately, the two scaling plots are of 
quite similar quality and with the system sizes considered here we cannot 
distinguish between the two critical scalings. 

4.2.2 Distribution of local susceptibilities 

An alternative analysis of the critical scaling in the disordered models is 
based on the study of the probability distribution function (pdf) of local 
linear and non-linear susceptibilities and how they behave when approaching 
the critical line. 

In the isolated random chain the pdf of local linear and non-linear sus- 
ceptibilities decay, for large values of the arguments, with a power law that 
decreases when approaching the quantum critical point. The inverse of this 
power is linked to the dynamic critical exponent and its divergence on the 
transition implies activated scaling. At a finite distance from the critical 
point the decay becomes sufficiently slow so as to lead to a divergent global 
susceptibility. One may wonder whether this phenomenon also occurs in the 
presence of dissipation. 

To try to give an answer to this question we started by studying the 
pdf of local linear susceptibilities of finite clusters {N < oo). In Fig. El we 
show p{\nxi) for = 4 and four values of the number of imaginary-time 
shces, A^^ = 32, 64, 128, 256, at a = 0.45 and AtT = 0.66. A finite value 
of Nt gives a finite bound to the maximum possible Xi- This figure can be 
compared to the left-most curve in Fig. 20 in the first reference in that 
has been obtained for = 4 and effectively with A^^ —>■ oo. This figure shows 
that the coupling to the external bath has not modified the form of the pdf 
in an important way. 

Finally, we analyzed the decay of the probability density of the local lin- 
ear susceptibilities when approaching the critical line. Each panel in Fig. 
displays these pdfs for four values of AtT and a fixed. Each set of data, 
represented with different symbols, were obtained for different A^ and choos- 
ing A^^ in such a way that the maximum in the Binder parameter is reached. 
The values of A^ are given in the captions. The straight lines are guides to 
the eye to indicate a possible power-law decay and the values of this power 
are given in the caption. 

Even if it is clear from the figures that there is a trend to a decreasing 
power law decay as the transition is approached, a functional relation between 
this power and the distance to the critical line is hard to establish beyond 
doubt from the pure numerical data. Analytic guidance, as was necessary 
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Figure 9: (Color on-line.) Probability distribution function of the local linear 
susceptibilities for a finite cluster with = 4 and (from left to right) four 
values of the number of imaginary-time slices, A^,- = 32, 64, 128, 256. a = 
0.45 and AtT = 0.66. 





Figure 10: Probability distribution of the single site linear susceptibility. 
AtT = 0.48 (top left), 0.43 (top right), 0.41 (bottom left), 0.39 (bottom 
right), a = 0.15 in all cases. The slope of the straight lines indicated in each 
panel are b = 4.55 (top left), 1.82 (top right), 1.20 (bottom left), 0.5 (bottom 
right). 
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to interpret correctly the numerical data for the isolated case, seems to be 
necessary here too. 

5 Conclusions 

With a careful analysis of the Monte Carlo data we succeeded in determining 
the phase boundary between the paramagnetic and ferromagnetic phases in 
finite and infinite rings of interacting quantum Ising spins in a transverse field 
and coupled to an external environment. While for finite number of spins 
there is no phase transition in the limit a — *• 0, in the thermodynamic limit 
one recovers the critical points of McCoy and Wu in the disordered case, and 
Onsager in the ordered problem. 

The method that we used allowed us to determine the Caldeira-Leggett 
critical value Oc = 1 for N = 1 with relatively little numerical effort. We 
used a similar analysis to find the phase transition of interacting finite clusters 
under the effect of an Ohmic bath. A similar analysis could be used to study 
the influence of other types of baths (sub-Ohmic, super-Ohmic) for which 
less analytical results are available. 

The coupling to the bath enhances the extent of the ordered phase, as 
found in mean-field spin-glasses jHl E] and in the ordered ferromagnetic Ising 
chain The Griffiths phase still exists as characterised by large distri- 
butions of the local linear and non-linear susceptibilities (the latter are not 
shown in the text). We expect to find similar results in higher finite dimen- 
sions and when frustration is included {i.e. when negative exchanges also 
exist). 

Our results for the type of scaling characterizing the phase transition 
for the disordered and dissipative infinite chain are not conclusive. Even if 
we have no evidence for the failure of the activated scaling of the isolated 
case when dissipation is included, we cannot exclude this possibility from 
our numerical data. It is important to notice that the activated scaling in 
the dissipation-less problem was found analytically and it was very hard to 
confirm numerically. It would be extremely interesting to extend Fisher's 
renormalization group approach to the case of the random chain coupled to 
an environment to give a definitive answer to this question. 

The systems we studied admit applications in a variety of fields. Quan- 
tum (spin) glass phases have been observed experimentally in several com- 
punds [32j; the role played by the coupling to the environment in their low- 
temperature anomalous properties has not been fully elucidated yet. Inter- 
acting spin systems are now becoming popular as possible toy models for 
quantum computers, with the Ising spins representing the qubits. It is ob- 
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viously very important in this context to establish the effect of the couphng 
to the bath. Another motivation for our study comes from the proposal in 
[T8] to describe non-Fermi liquid behavior in certain /-electron systems with 
a similar model (see also JHI)- Finally, while classical phase transitions are 
well understood by now, the same does not apply to quantum phase transi- 
tions at zero temperature jSl ESI • ■^o'^ these reason we believe that this 
(and related ones) are interesting problems that deserve further study. 
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